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Abstract 

We present a calculation of the Vector Form Factor at the next-to-leading order in the 
I/Nq expansion, within the framework of Resonance Chiral Theory. The calculation 
is performed in the chiral limit, and with two dynamical quark flavours. The ultravi- 
olet behaviour of quantum loops involving virtual resonance propagators is analyzed, 
together with the kind of counterterms needed in the renormalization procedure. Using 
the lowest-order equations of motion, we show that only a few combinations of local 
couplings appear in the final result. The low-energy limit of our calculation reproduces 
the standard Chiral Perturbation Theory formula, allowing us to determine the reso- 
nance contribution to the chiral low-energy couplings, at the next-to- leading order in 
1/Nc, keeping a full control of their renormalization scale dependence. 



1 Introduction 



We have at present an overwhelming experimental and theoretical evidence that the SU{3)c 
gauge theory correctly describes the hadronic world [1]. This makes QCD the established 
theory of the strong interactions. Nevertheless, its non-perturbative nature at long distances 
is still challenging our theoretical capabilities. Since the hadronization procedure is not 
understood, an effective field theory description [2] in terms of hadronic degrees of freedom 
is required in the low-energy regime. 

Below the heavy quark thresholds, QCD is properly described considering only the light 
quarks, with masses much smaller than the dynamical QCD scale Aqcd- One can then study 
the massless QCD case and consider the mass term as a perturbation. With Uf massless 
quarks, QCD has a chiral SU{nf)ji ® SU^nf)^ symmetry which is spontaneously broken 
to SU{nf)L+R- This generates {n'j — 1) Goldstone bosons, which can be identified with the 
multiplet of light pseudoscalars; their small masses being proportional to the explicit breaking 
of chiral symmetry generated by the quark masses. The Goldstone nature of the pseudoscalar 
bosons implies strong constraints on their interactions, which can be most easily analyzed 
on the basis of an effective Lagrangian organized as an expansion in powers of momenta 
(derivatives) and quark masses over the chiral symmetry breaking scale ~ 471^,^ ~ 1.2 
GeV. The resulting Goldstone effective field theory, Chiral Perturbation Theory (xPT) [3-5] , 
has achieved a remarkable success describing the low-energy dynamics of QCD [6,7]. 

In the resonance region one must introduce a different effective field theory with explicit 
massive fields to describe the degrees of freedom associated with the mesonic resonances [8,9]. 
Although chiral symmetry still provides stringent dynamical constraints, the usual xPT 
power counting breaks down in the presence of higher energy scales. Therefore, one needs 
another expansion parameter to organize the effective Lagrangian. The limit of an infinite 
number of quark colours [10] turns out to be a very useful tool to understand many features 
of QCD and provides an alternative power counting to describe the meson interactions [11]. 
Taking Nc — > oo, with agNc fixed, there exists a systematic expansion of the SU{Nc) gauge 
theory in powers of 1/Nc, which for Nc = 3 provides a good quantitative approximation 
scheme to the hadronic world [12]. Assuming confinement, the strong dynamics at Nq — > oo 
is given by tree diagrams with infinite sums of hadron exchanges, which correspond to the 
tree approximation to some local effective Lagrangian. Hadronic loops generate corrections 
suppressed by factors of l/Nc- 

The short-distance properties of the underlying QCD dynamics impose strong constraints 
on the couplings of the hadronic effective theory [9,11]. The infinite sums of meson exchanges 
contributing to any given Green function should obey the right QCD behaviour G(t) ~ 
at large momenta. This requirement excludes resonance interactions with large number of 
derivatives, explaining the phenomenological success of the usual lowest-order (in derivatives) 
approximations. Moreover, it implies stringent correlations among those resonance couplings 
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associated with the highest powers of momenta. The usual xPT expansion is of course 
recovered at very low energies, when the resonance Green functions are expanded in powers 
of momenta over the resonance mass scale. 

In spite of its many dynamical simplifications, QCD at Nq — > oo is still a very involved 
theory and some approximations arc called for. Usually one cuts the infinite tower of reso- 
nance exchanges to a finite number, taking only into account those meson states which are 
relevant at the physical energy scale. This is meaningful since the contributions from higher- 
mass states are suppressed by their corresponding propagators. However, it introduces back 
a momentum expansion regulated by inverse powers of the heavier resonance masses which 
have been integrated out. The problem can be formally avoided taking the limit Mr oo 
for all resonance states not included in the effective theory. This gives a well-defined ap- 
proximation with a clear physical meaning: one is assuming that the QCD short- distance 
operator product expansion provides an acceptable description at energies above the last in- 
cluded mesonic state. The imposed short-distance constraints are nothing else than matching 
conditions between the low-energy effective field theory and the underlying QCD dynamics. 

The most drastic and simplest scheme is the so-called Single Resonance Approximation, 
which only considers the contributions from the lightest meson with any given quantum 
numbers [8,9,13,14]. The short-distance QCD constraints determine in this case all hadronic 
parameters in terms of the pion decay constant F.,^ and the two masses of the vector and 
scalar multiplets. My and Ms [11]. This gives a very successful description at energies below 
the scale of the second resonance multiplets. Since there is an infinite number of Green func- 
tions, it is clearly not possible to satisfy all matching conditions within the single resonance 
approximation. A useful generalization is the Minimal Hadronic Ansatz, which keeps the 
minimum number of resonances compatible with all known short-distance constraints for the 
problem at hand [15]. The resonance contributions to some 0{p^) xPT couplings have been 
already analyzed in this way, by studying an appropriate set of three-point functions [16-18]. 

The Nc oo limit of the resonance chiral theory (R^T) has been investigated in many 
works [8,9,11,13-18] and a very successful leading order phenomenology already exists [11]. 
However, we are still lacking a systematic procedure to incorporate next-to-leading contribu- 
tions in the 1/Nc counting. Up to now, the effort has concentrated in pinning down the most 
relevant subleading effects, such as the resonance widths which regulate the corresponding 
poles in the meson propagators [19-21], or the role of final state interactions in the physical 
amplitudes [19-27]. These unitarity corrections are generated through Goldstone loops and 
therefore are suppressed by l/Nc factors; nevertheless they may be strongly enhanced by 
large infrared logarithms. The combined constraints of analyticity and unitarity make pos- 
sible to perform appropriate resummations of chiral logarithms, which describe the leading 
1/Nc corrections in the resonance region. 

Quantum loops including virtual resonance propagators are a major technical challenge 
which has never been properly addressed. A first step in this direction was the study of reso- 
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nance loop contributions to the running of the xPT couphng Lio(/i), performed in Ref. [28], 
which however didn't attempt any analysis of the induced ultraviolet divergences and their 
corresponding renormalization.^ Clearly, at the one-loop level the massive states present in 
RxT generate all kind of ultraviolet problems which are not yet understood. A naive chiral 
power counting indicates that the renormalization procedure will require higher dimensional 
counterterms, which presumably could generate a problematic behaviour at large momenta. 
Therefore, it will be necessary to perform a careful investigation of the constraints implied 
by the short-distance properties of QCD at the next-to-leading order in 1/Nc- 

A formal renormalization of R^T at the one-loop level appears to be a very involved task, 
which requires the prior analysis of several technical ingredients [34]. In order to gain some 
understanding on the ultraviolet behaviour, it seems worth to perform first some explicit one- 
loop calculations of well chosen physical amplitudes. In the following, we present a detailed 
investigation of the pion vector form factor (VFF) at the next-to-leading order (NLO) in the 
1/Nc expansion. This observable is defined through the two-Goldstone matrix element of 
the vector current: 

{n+{p,)n-{p,)\^-{uru-drd) |0) = ^(g^) (pi - , (1) 

where q'^ = {pi + P2)'^- At very low energies, the VFF J^{q^) has been studied within the 
xPT framework up to C(p^) [35-37]. R^T and the 1/Nc expansion have been also used 
to determine J^{q^) at the p meson peak, including appropriate resummations of subleading 
infrared logarithms [19-22]. 

We will simplify the calculation working in the two flavour theory and taking the massless 
quark limit. Therefore, we will assume a chiral U{2)l ® U{2)r symmetry group. The small 
effects induced by the U{1)a anomaly will be neglected, because they are not going to be 
relevant in our discussion. The isosinglet pseudoscalar can only appear within loops, and 
the numerical correction generated by its non-zero mass could be taken into account in a 
straightforward way, together with the finite quark mass effects which we are ignoring. 

In the next section we will briefly describe the R^T Lagrangian. We will adopt the 
single resonance approximation and will only consider the minimal set of resonance couplings 
(linear in the resonance fields) introduced in ref. [8], supplemented with those counterterms 
required by the renormalization procedure. The renormalization of the relevant one-particle- 
irreducible Feynman diagrams will be discussed in section El and the final results of our 
calculation will be collected in section |3 Sections El and IHl analyze the behaviour of the 
computed vector form factor at low and high energies, respectively. We will finally summarize 
our findings in section [71 Some technical details have been relegated to the appendices. 

^Quantum loops involving massive states have been only analyzed within explicit models with additional 
symmetries. For instance, the gauge structure advocated in the so called "Hidden Local Symmetry" descrip- 
tion of vector resonances [29] implies a much simpler ultraviolet behaviour [30] . Chiral loop corrections to 
some vector resonance parameters have been also studied [31,32] within the context of "Heavy Vector Meson 
xPT" [33], which adopts the Mj^ oo limit to guarantee a good chiral power counting. 
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2 The Lagrangian of Resonance Chiral Theory 



We are going to work within a U{2)l ® U{2)r chiral theory, containing a multiplet of 4 
pseudoscalar Goldstone bosons, 

parameterized through the unitary matrix = exp |i $/(-\/2 The Goldstones couple 
to massive U{2) multiplets of the type V(l—), S'(0++) and P(0-+), with a field 

content analogous to the one indicated in 

Our starting point is the R^T Lagrangian C{u,V, A, S, P) introduced in Ref. [8]. It 
contains the 0{p^) xPT Lagrangian [4,5], 

C2>c = -^{uX + (3) 
the kinetic resonance Lagrangians, 

C2z{R = V,A) = -\ly^R^^S/,R-^^-\MlR^,Rn. 

C2ziR = S,P) = \{\/>'R\/^R-MlR^), 
and 0{p'^) interactions linear in the resonance fields: 



(4) 



(5) 



aa[a(i++)] = ^{A,jn, 

As[5(0++)] = CdiSu^u'^) + c^{Sx+), 

C2P[P{0-+)] = td„^{Px-). 

The brackets (...) denote a trace of the corresponding flavour matrices, and the different chi- 
ral tensors follow the notation defined in Ref. [8]. They include external vector, axial- vector, 
scalar and pseudoscalar sources {v^, a^, s, p) to generate the corresponding Green functions. 
Following this reference, we describe the vector and axial-vector resonances with the anti- 
symmetric field formalism. In the chiral limit and neglecting external scalar or pseudoscalar 
sources x± = 0. 

The Lagrangian C{u, V, A, S, P) obeys the correct Nc counting rules. The different fields 
and the masses and momenta are all of them whereas all couplings (F, Fy, Gy, Fa, 
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Cd, Cm and dm) are of O i^^/Ncj- In the limit Nc oo, one can determine all parameters 
in terms of F, My and Ms [11]. The short- distance QCD behaviour of the vector, axial- 
vector [9] and scalar [24] form factors, together with the constraints implied [9, 11, 14] by 
the superconvergence properties of the vv — aa [38] and ss — pp two-point functions at large 
momenta, imply the relations [11]: 

^ = V2Gv = FA = 2cm = 2cd = 2V2dm = F, (6) 

Ma = V2Mv , Mp = V2Ms (1 - Sf^^ . (7) 

The last identity involves a small correction 5 ~ SirasF'^/Mg ~ 0.08 a^, which can be 
neglected together with the tiny effects from light quark masses. 



2.1 Subleading Lagrangian 

The one loop calculation of the VFF with the previous Lagrangian generates ultraviolet 
divergences which require counterterms with a higher number of derivatives. We will only 
include the minimal set of chiral structures needed to renormalize the VFF calculation. We 
expect their corresponding couplings to be subleading in the 1/Nc expansion, since they are 
associated with quantum loop corrections. 

We need to include the following 0{p^) and 0{p^) Goldstone interactions: 

Ax = ^ ( /r K, nj\ ) - I,, ( V^u.V^u, ) , (8) 

C^X = ^ C51 ( V^/r [/^MP> U,]) +1 C53 ( V^/r [/i.p, U"] ) . (9) 

Chiral invariance forces these terms to have structures contained in the corresponding xPT 
Lagrangians [4,39]. We use a tilde to denote the R^T couplings in © and 0, which are 
different to the ones with the same names (without tilde) in xPT. For instance, the xPT 
coupling ^6 (-^9 in the three flavour case) is dominated by a contribution from vector-meson 
exchange and is of 0{Nc), while the corresponding R^T coupling does not contain this 
contribution and is of 0{1). 

The 0{p'^) term containing the structure V^u^ does not contribute to the tree-level 
calculation of the VFF; nevertheless, it is needed to renormalize the Goldstone self-energies. 
At 0{p^), only the combination of couplings = 4F^ (C53 — C51) is going to be relevant for 
the VFF [39]. Including the Lagrangians (jHI) and Q, the tree level calculation of the VFF 
gives the result: 

^(9 ) = 1 + -p^w— ^ - ^^Y^ + ^viY^- (10) 
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Figure 1: Leading order contributions to the VFF. 



The requirement that the VFF should vanish at gr^ — > oo imphes the following conditions at 
leading order in l/Nc'- 



FvGv ^ , 4 = , ry^ = AF' (^53 - ^51) = 0. (11) 

Therefore, the couplings Iq/F'^ and Vy^jF'^ are of subleading order in the 1/A^c expansion, 
i.e. C(l/iVc), as expected on pure dimensional grounds. 

The renormalization of Green functions including resonance fields forces the presence of 
the following additional counterterms: 



Uz= ^ ( v^y^^ {v., vn V ) + ^ ( {v., V J v^- {v^ V"} 



(12) 



+ ^ ( {v^ V"} v^^ {v., V J V,, ) , 

= Xp^ { i;,VVr ) + XfA V^, {W, V J f^" ) , (13) 
Ag ^iXcA {V^ V^} V-" [u,, Uo]) +t Xg, { V'^'' K^, K] ). (14) 

The quadratic Lagrangian C^z is needed to renormalize the vector self-energy. Actually, 
only the sum of couplings Xz = Xz^ + Xz2 + Xz^ is relevant for this purpose. The 

renormalization of the v-V two-point function involves the sum of /^4f couplings Xp = 
Xp^ + Xf2- Finally, the three-point function with one external vector resonance and two 
Goldstone legs is renormalized by £4(5 through the combination Xq = Xg2 — \Xg^- The 
dimensions of the couplings are [Xz] = E'"^ and [Xp] = [Xq] = E~^. 

At the NLO in I/Nq, these countcrtcrm Lagrangians only contribute through tree-level 
diagrams. One can then use the leading order equation of motion (EOM) 

V^'VpVP" -VVpYP^" = -M^V^"" - - ^ K,-""] (15) 
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to reduce the number of relevant operators. The Lagrangians fll2|) . ()13p and ()14|) take then 
the equivalent forms: 



^EOM ^ XzMy ^ y^^y^^ ^ ^ ^^^^ ( ^../r ) + ^^^^^ ( V,^ [u^ ) 



(16) 



= -x,Mi { v;./r ) - ^^S^ ( /r «.]) + ■■•, (17) 

= -2^ XgM2 ( V^"- n,] ) - z X^Fv. ( /f [m^, «.]) + ■••, (18) 

where the dots denote other terms which are not relevant for the VFF calculation at this order. 
The derivatives acting on the vector resonance fields have been traded by the heavy mass 
scale My and/or derivatives acting on the Goldstone fields, giving rise to the usual tensor 
structures of the xPT Lagrangian. Therefore, the effect of the counterterm Lagrangians C^z, 
and CiG is just equivalent to the following shift in the couplings at the next-to- leading 
order in 1/Nc'- 

If = 4 + 2XzFvGv - 2y/2 XpGv - ^V^XcFy, 
Ff = Fv + 2XzM^Fv - 2^2 XpM^, 

Gf = Gv + 2XzM^Gv - AV2XgM^, (19) 



V ; 



-eff ~ 

V2 — ''^V2 



Thus, since if ^ ^ {Ml'f^ ^ r-. 0{l) and Ff ~ Fy ~ Gf ~ Gy ~ 0{^/N^), a 
consistent 1/Nc counting requires that Xq and Xp are of 0{1/ ^/ Nc) and Xz of 0{1/Nc)- 



3 Renormalization of Quantum Loops 

The renormalization procedure follows very systematic and precise steps in any well defined 
field theory. First of all, the two-point Green functions must be renormalized. Later the 
three-point Green functions and so on. For the VFF up to NLO in 1/Nc only the two- and 
three-point Green functions will contribute. The corresponding renormalizations for the one 
particle irreducible diagrams (IPI) at one loop are given in the next subsections. 
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We will adopt the MS — 1 scheme, usually employed in xPT calculations, where one 
subtracts the divergent constant 

Aoo = ^/^+7i?-ln47r-l. (20) 
a — 4 

However, we will impose the on-shell condition to renormalize the pion self-energy. This 
simplifies the calculation of physical amplitudes with external pions. Since we work in the 
massless quark limit, the Goldstone tadpoles will not give any contribution. The precise 
definition of the relevant Feynman integrals with one, two and three propagators and some 
useful antisymmetric formalism technology are relegated to appendices A and B. 



3.1 Pion self-energy 




Figure 2: One-loop diagrams and local contributions to the pion self-energy. 



The diagrams contributing to the pion propagator are shown in Fig. |21 The kinetic 
Lagrangians C2z{R) in Eq- © generate additional tadpole topologies with one resonance 
propagator, but they are identically zero even with massive pions. The divergences of 0{p^) 
are reabsorbed through the wave-function renormalization vr^^^ = (1 -|- 5Z^)2 vr*^'''', being vr^'^'' 
and TT*^^) the bare and renormalized pion fields respectively. In the on-shell scheme. 



1G\ 



?>Ml 



Aoo + In ■ 



/i 
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+ 



4d 



Ml 



Aoo + In - - 



There are also divergences of O(p^) which renormalize one of the couplings in £4^: 
The renormalized pion self-energy takes the form 



Gy + 2 Aoo 
P 32^ 



p 



167r2F2 



2Gl r 



In- 



2 ^ 



P 

Ml, 



(21) 



(22) 



+ 



4d 



In^ + 



P 

Ml 



(23) 



where the function 

4.(.) = (i - ' 



X 



X 



1 - - In 1 - X - 1 + - 



X 



^0 (n + 2)(n + 3) 



(24) 



contains finite and scale-independent contributions. 



3.2 Rho self-energy 




Figure 3: Rho self-energy. 



The one loop p self-energy contains only an 0{p'^) divergence, which renormalizes the 
coupling Xz of the NLO resonance Lagrangian: 



Xz = X'M+5Xz{p) 



SXzifi) 



2Gl 



Xr 



2 772 



Thus, the vector mass and wave-function are not renormalized: 

5Ml = ; 5Zy = 0. 

The renormalized p self-energy then becomes: 

where the antisymmetric tensor structure A^^'''"{q) is defined in appendix B and 



2G 



V 



F' F' 



-Bo(q^/p^) + -^-^ 



(25) 
(26) 
(27) 

(28) 



3.3 {v^" VP") IPI vertex 

The IPI amputated diagrams (at NLO) connecting an external vector quark current to an 
outgoing vector resonance are shown in Fig. |3J The one-loop contribution brings an 0{p'^) 
divergence which gets reabsorbed through the following renormalization of the coupling Xp: 

Xp ^ X^M + 6Xp{p) ■ 6Xp{p) = . (29) 
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6Fy,5Zy,Xp 

(go V 



Figure 4: Diagrams contributing to the {v^V^") Green function at NLO in l/N^ 



c- 



Since there are no divergences of O(p^), the lowest-order couphng Fy remains unchanged: 



5Fy = 0. 

The renormahzed vertex function takes the form 



1 - 1 
-Bo(q'^/a^) + ^ 



(30) 



(31) 



where the first term is the leading order contribution. The massless two-point function 
Bo{q^ / /i^) is defined in appendix A and the antisymmetric tensor structure X^^ in appendix B. 



3.4 {V^^TiTi) IPI vertex 




Figure 5: NLO diagrams contributing to the three-point Green function V^'^ — > vrvr. 

The IPI amputated diagrams connecting a vector resonance with two outgoing pions at 
NLO in 1/Nc are shown in Fig.|3| The loop diagrams generate 0{p'^) and 0{p'^) divergences, 
which renormalize the couplings Gy and X^, respectively: 

Gv = Gy{fi) +SGv{n) ; SGvifi) = Gy 

(32) 



3Mt 



2Gl 



b p2 



levr^F^ 
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<5XG(/i) 



2^ 



-^ + ^-2 



Xr 



(33) 



The wave-function renormalization of the external vector and pion legs amounts to a global 
factor (6Zt, + ^6Zv^ multiplying the lowest-order contribution {dZy = at NLO). Taking 
this into account, one finally gets the finite vertex function 



^r;^(pi,P2) 

where 

A^(g^/x2) = 



F2 



+ 



2G? 



6 



V 



In- 



Mi 



2Gl (Ml 



+ 5 



Ml 



In 



MlA^fMl 



+ 



V 



647r^ 



3^ - 1 



3M|4c| 

"T " • 2 77^2 "1" 



+ 



^Co(g^O,0,M^, 



647r^ 



2887r' 



2Gl 



4d 



i?2 + ^2 



- 2 



+ ^Co(g^o,o,M|) 



+ 



g2 + 2 



m: 



(35) 



The three-propagator integral Co(g^, , M^^, M^) is defined in appendix A. 



3.5 (viiTttt) IPI vertex 



The divergences generated by the IPI loop diagrams shown in Fig. get reabsorbed through 
the renormalization of the pion wave function 6Zt^ and the 0(p'*) and 0{p^) couplings and 



'V2- 



4 ^ + 6l,{^^) 
rv2 = rv^ifi) + 6ry2{n) 



54(/i) = 



F'K 
9Qr 



1 1 

— + — 



M^ 



M', 



The resulting finite correction to the lowest-order pion form factor, 

A^(g2)iPi = AJ-^=* + Ajj2 + A^J + AJ-^ + + A J", 
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■p 



(36) 
(37) 

(38) 




X 



Figure 6: IPI diagrams connecting an external vector current and two outgoing pious. 



contains contributions from tree-level counterterms, 



2Gl 



Ml 



Ml 



4cj 



Ml 



3 In + - I + zr^ 

,4 
_p4 ' 



2 



and loop diagrams with internal Goldstone bosons (first diagram in Fig. | 



|,{i5o(gV/^^) + 



1447r^ 



and vector. 



(39) 



(40) 



V 



2Gi 1 



F 



-2 -Co(g^O,0,M^ 



2 ^ 



5M^ 



+ Co(g^M2,M2,0) 



g 2 



i?o(g7/) 



^ + 2M2 + 



12 



5o(g',M^ 



2M2 + ^ 



6 6M, 



+ 



2G^ 



V 



+ 



2M, 



5g 
l2 



2\ n 



(41) 



+ 



IGtt^F^ 



-1 + 



12M^^ 



2M^ 



axial-vector, 



5o(g^Ml 



2g^ ^ 



+ 



2Gl 



'Ml 



+ 1 - 



19g^ ' 



2M2 + ^- ^ 



+ 



'A 

^4 



ml 



+ 



Ml ,3 



levr^F 
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2Ml 



6Ml 



3271^^2 



727r2F2Ml ' 



(42) 



12 



scalar, 



^^{-C^o(g^O,0,M|) 



M| Ml 

2 I 

g 2 



+ Co(g^M|,M|,0) 



M| 







M| 1 




12 


+ 16 A' J 



2887r^F2 



1-1 ^ 



B^{q\M. 



'2Ml 



4c^ 



'Ml 



Ml 



3 12, 



(43) 



Ml 



4c^ 



and pseudoscalar resonances, 



2Ml ^ 



12M|^ 



+ 



6M|j 



6 
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In 



M: 



(44) 



4 Vector Form Factor 




The basic topologies contributing to the VFF are shown in Fig.|3 in terms of the NLO IPI 
diagrams computed in the previous section. The internal p line denotes the dressed vector 
propagator, including the self-energy correction ()28|) which regulates the p pole. Taking this 
self-energy into account, the LO contribution takes the form: 



ILO 



FvGl,{p) 



M^-q^- S(;(g^ 



(45) 



The topology in Fig. [7|a generates the following NLO correction: 



M^-q'- S(;(g2 




6 



^o(gV/i' 



1447r' 



-2F^XU^^) 



(46) 
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The second topology (Fig. [7|b) brings the contribution: 



(47) 



where Ar(g^,yU^) is given in Eq. (j35|) . Finally, Fig. [3c denotes the IPI correction AjF(g^)ipi 
in Eq. Adding all contributions together, one gets the VFF at NLO: 

^(g2) = ^(g2)^Q + A^(g2)^ + A^(g2)G + A^(g2)ipi. (48) 



Using the large- A^^c relations © in the NLO terms, the result can be written in the form: 



M2 -g2_S^(g2 



(49) 



where 



The constants 



Bi 



6v = 



l + 5v + 2M^X ~ AT{q^), 
Qiq^) -6v-2{M^ + q^) X . 



(50) 



X 



-i-Ar(o,/x2), 



[x];(/i) + 4x^(/i)] 



and the functions Sy(g2), 



and 



Af(g2) ^ Ar(g2,/x2) - Ar(0,/x2) 



(51) 
(52) 

g{q') ^ A^(g2)ipi + Af(g2) ^ G(g^/i2) - Ar(0,/i2) (53) 

are independent of the renormalization scale fi. The subleading R^T couplings Xp[ij,) and 
Xq[ix) only appear through the constant X, while X^(/i) is also present in the function 
S^(g2). At g2 = 0, Af (0) = ^(0) = S(;(0) = 0. Therefore J^(0) = 1 as it should. 

Some IPI diagrams (Figs. IHla and IHle and the V terms in Figs. IHlb and Elc) have 
a corresponding reducible counterpart involving a vector propagator. The combination of 
both contributions can be then incorporated in A{q'^). The function (^(g^,/!^) contains the 
corrections generated by the other IPI diagrams (Figs. Eld, Elf, the S term in Fig. Elb, 
the S, A and P terms in Fig. Elc and the and pieces in Fig. Elg)- Subtracting their 
contribution at q^ = 0, which contains the dependence on the renormalization scale fi, 



G{0,f,') = Ar(0,/i2) 



2 



2 +4 



/i 2 



(54) 
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one gets: 



+ F2 



Ml 
2 

1-2 



+ 



Co(g^M|,M|,0) 



Ml 



_s_ 

,2 



q' 3 + 4 + _ 



+ 



Bo{q\M\ 



Ml' 



Ml 



+ 



+ 



Bo{q\M\) 



-2Ml - — + ^ 

^ 6 6Ml_ 



+ 



Bo{q\Ml) 



2Ml 



+ 



'L 
6 



1 



+ M| 3 , . 1 . 



6 



1 ln^ + -ln- 

12 + o ,.2 



Mi 



In 



1 , Ml 1 , 1 , M|, 
- In In In + 

2 12 /i^ 6 /i^ ^ 

Ml A ( I 1 



3 12 



(55) 



^-167r2£^(/.) 



/i^ +3 l^M^ ^M^ 



+ 



^? M] 



+ 



967r2 



5 Low-Energy Limit 



At very low energies, <^ M^, the resonance fields can be integrated out from the effective 
theory. One recovers then [8] the standard xPT Lagrangian, which leads to the following 
resuh for the VFF [35,36]: 



•^xPT(g^) = 1 - |. {W + ^ 



In 



/i 



(56) 



+ 



V2 



967r' 



In -- 



(2£[ - + 40 (/i) + O I 



AT, 



The Taylor expansion in powers of q^ of the R^T prediction (j49p reproduces the xPT formula 
(jSni), as it should. The coefficient of the O [q^ In (— g^//^^)] term satisfies the known \axge-Nc 
equality [8,11] (2£i - £2 + 4) = i^^(l - 5M|/M^)/(2M|). The non-logarithmic C»(gO and 
0{q^) terms relate the low-energy chiral couplings 1^ and with their R^T counterparts 
^6 and ry2. 



M^ 



967r^ 



In 



9 + 3 In r. 



13 



15 



+ w + - 



Mt 



4 , M, 
- m 
3 /i 



1 , M 



f --In 

2 2 



7^ 



Ml , 

^ In^ 

Ml ^? 



11 M| 
^ 18 ^ 2M^ 



(57) 



y2 



+ 



967r^ 




2M^ 



2M|^ 



In 



Mi 



Ml 



In- 



167 



17 



lOMl 



3M| 



+ 



17 



20M| 



+ 



10M|, 



(58) 



Notice that the combination of subleading R^T couphngs X does not appear at 0{p^). 
Therefore, the relation (|57p adopts the same form in terms of the effective couphngs defined 
in (Uni), i.e. if'''{^^) - Gf'''ii^)/{Miy''''{i^) = £^(/i) - Fy G{;(/i)/M2 . As shown in (|SHI), 
this is no longer true at 0(p®); nevertheless, the explicit dependence on X — X'^^fi) present 
in ry2if^) can be reabsorbed into the leading term, through the use of the effective couplings. 



i.e. 2(/i) = F'Ff '^(/x)/(M^)^ff''^(/i) + rf/ + 



Eqs. ()57p and (j3H|) contain the well known lowest-order predictions for the two xPT 
couplings: = —Myry2/F'^ = —F'^/My. Moreover, they give their dependence on the 
renormalization scale at the NLO. The running of the renormalized couplings [^e (a*)' '"y2(A'')] 
and [ifi{fJ'),ry2if^)] is different, because their corresponding effective theories have a very 
different particle content. 



The fi dependence of a given coupling "g" can be characterized through the logarithmic 
derivative 

— = (-59) 
From Eqs. ()36|1 and ()37|1 one gets the running of the R^T couplings: 



7^ 



2 

3 




^ Ml 



2Mt 



(60) 



Eqs. (j57j) and (j58p give then the dependence on the renormalization scale of the corresponding 
XPT couplings: 

1 F2 / 5 1 \ 

X. = -o ' X = — I T7T - 772 I • (61) 



Ml 



These values are in perfect agreement with the low-energy results of refs. [4,5,36,39]. The 
running of the 0{p^) coupling ry2{^)/F^ receives of course additional 2-loop contributions 
which are oiO{l/Nl). 
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The rigorous control of the renormahzation scale dependences allows us to investigate the 
successful resonance saturation approximation [8] at the NLO. The xPT couplings and 
have been phenomenologically extracted from a fit to the VFF data at low momenta. 
This determines [36] the scale-invariant combination 



(/i) - log^ = 16.0 ±0.5 ±0.7, (62) 
and 

r^2(M^) = (1.6 ±0.5) • 10-^ (63) 

Inserting these numbers in Eqs. (j57jl and (f58|l . one can estimate the corresponding scale- 
invariant combinations of NLO couplings in R^T: 

4-?;M-5^1og^-j^4.. (64) 



F 



_ 2F' 



rv2 ^ rl.,{f,) ± V74 (^v + 2M^ [x - (/.)]) - '7 log , (65) 

My ^ L OZTT fi 

where 7^^ = -1/(6^^). Taking F = 92 A MeV, My = 770 MeV and M5 = 1 GeV, one 

gets 4 = (-0.2 ± 0.9) ■ 10"3 and = (-0.2 ± 0.5) ■ 10-^ while a larger value of the 
scalar resonance mass Ms = 1.4 GeV shifts the 0{p^) coupling to Iq = (—0.9 ± 0.9) ■ 10~^, 
without affecting at the quoted level of accuracy. These numbers should be compared 
with the laxge-Nc predictions for the xPT couplings iel^^^oo = —F'^/My = —0.014 and 
rygliVc^oo = F^/My = 2.1 ■ 10"'^. Put in a different way, the hypothesis = fy^ = 
generates excellent predictions for (yu) and ry2(/^) at any scale fi. 



6 Behaviour at Large Energies 



At large momentum transfer, the relevant renormahzation scale invariant functions take the 
forms: 



1 



16n^F 



2 772 



6 



1 1 

± 



Ml 



In 



^ln^±y-167r^4^(^) 



+ on 



2 ^V2if^) 



Af(g2) 



V 



Wn^F 



2 Tp2 





In 2 


-lln^ 






2 



9 m: 



± 



s 



IT 



6 ■ 4 ■ 4M^ 



,2 ' 



(66) 
(67) 
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The p propagator makes the A{q^) piece of the VFF well behaved when oo. However, 
the IPI contributions generate a wrong behaviour Q{q^) ~ In {—q^/li'^) in the B{q'^) term, 
which cannot be eliminated with a local contribution. The problem originates in the two- 
resonance cut which has an unphysical growing with momenta. Although our leading R^T 
Lagrangian (0) only incorporates couplings linear in the resonance fields, the kinetic resonance 
Lagrangian (@]) introduces some bilinear interactions through the chiral connection included 
in the covariant derivatives. Their couplings are fixed by chiral symmetry and give rise to the 
IPI diagrams in Figs. IHlc, IHld andlHlf- Obviously, these are not the only interactions bilinear 
in the resonance fields even at large- A^^c [17,18,40,41]. Therefore, it is not surprising that our 
calculation is unable to find the correct behaviour at large energies for those contributions 
with two intermediate resonances. 

The contributions with an internal vector propagator in diagramsElb andlHlc give us some 
hint about which pieces could be missing in our calculation. These two diagrams combine 
with a reducible contribution of the type [Zlb: the IPI (V^i/VrTr) vertex in Fig. Elb. The 
three contributions contain identical loop functions and their sum generates a global factor 
My /{My — g^), which suppresses the large-g^ behaviour. Thus, these corrections have been 
included in the term A(g^). 

It seems natural to conjecture that the remaining IPI contributions with two-resonance 
cuts should combine with the corresponding reducible topologies, including {VRR) and 
{v^RR) vertices, to generate the final propagator suppression: 

Gte^) - I^TTf^Gtf). (69) 

The needed Lagrangian takes the form 

(-"vRR = ^ A^^^ ( V^^ V^S ) + z A^^^ ( V'" V^P V.P ) + ■ ■ ■ (70) 

Our conjecture fixes the new chiral couplings in the large-Nc limit. It would be interesting 
to analyze the contributions of this Lagrangian to appropriate Green functions, following 
the work of refs. [16-18], and check whether the couplings predicted by the corresponding 
short-distance QCD corrections agree with our naive conjecture. In appendix C, we show 
two simple examples where the presence of the propagator suppression can be demonstrated 
in a rather straightforward way. 

The behaviour at large energies is also constrained by unitarity requirements. Moreover, 
the local contributions can be forced to vanish at large q^ by taking appropriate values of the 
RxT couplings. Probably, this could allow us to determine the scale invariant constants 
and fy2. We plan to investigate all these points in forthcoming works. 
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7 Summary 



The one-loop analysis of the VFF has shown a series of interesting features. As expected, loop 
diagrams with massive resonance states in the internal lines generate ultraviolet divergences, 
which require additional higher- dimensional counterterms in the R^T Lagrangian. Since 
these counterterms give rise to tree-level contributions which grow too fast at large momenta, 
their corresponding couplings should be zero at leading order in the large-A'^c' expansion. 
Thus, one can establish a well defined counting in powers of I/Nq to organize the calculation. 

The formal renormalization is completely straightforward at one loop. One can easily 
determine the dependence of all relevant renormalized couplings. Moreover, the final result 
is only sensitive to some combinations of the chiral couplings. In fact, using the lowest-order 
equations of motion, one can eliminate most of the higher-order couplings. Their effects get 
then reabsorbed into redefinitions of the lowest-order parameters. 

Expanding the result in powers of g^/M^, one recovers the usual xPT expression at low 
momenta. This relates the low-energy chiral couplings and with their corresponding 
RxT counterparts and f^g- The rigorous control of the renormalization scale dependences 
has allowed us to investigate the successful resonance saturation approximation at the next- 
to-leading order in I/Nq. The assumption = = generates excellent predictions for 
£g (/x) and ry2{lj) at any scale ii. 

At high energies, we have identified a problematic behaviour which originates in the 
two-resonance cuts: they generate an unphysical increase of the VFF at large values of 
momentun transfer. This is not surprising, since there are additional contributions generated 
by interaction terms with several resonances, which have not been included in the minimal 
RxT Lagrangian. These new chiral structures should be taken into account to achieve a 
physical description of the VFF above the two-resonance thresholds. The short-distance 
QCD constraints can be used to determine their corresponding couplings. 

Our calculation represents a first step towards a systematic procedure to evaluate next-to- 
leading order contributions in the 1/Nc counting. More work in this direction is in progress. 
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Appendix A: Feynman Integrals 

The calculation involves the following Feynman Integrals: 



Aoo + In ■ 



M2 
7 



(A.l) 



{2Try {k' + ie- M^) [{q - k)' + ie- M^] 



(A.2) 



and the finite function 

C,{q\MlMlMl) ^ 



In—^ 



Ml 



Ml - Ml I,' 



+ J{q',Ml,Ml), 



(A.3) 



2\ ' 



t{27iy [{pi - ky + te - Ml] [(p2 + ky + te - Ml] {k' + te - Ml) 

with q = Pi + P2 and, with massless outgoing pions, — P2 — 0. The divergences are 
collected in the factor 

,d-4 



A. 



2 11° 



d-A 



+ 7b - In 47r - 1 , 



(A.4) 



being 7^ ~ 0.5772 . . . the Euler constant and ji the renormalization scale. 
The two-propagator integral contains the finite function 



J{q\MlMl) 



1 



327r2 



2 + 



Ml - Ml Ml + Ml 



Ml - Ml 



ln4 



(A.5) 



\'l\q\ Ml Ml) ,^ / [g2 + X^^q^ Ml Ml)f - {Ml - M, 



In 



ir 



r \ [q^ - AV2(g2, Ml Ml)f - {Ml - M^ 

with X{x, y, z) = x"^ + y'^ + z"^ — 2xy — 2xz — 2yz . Some useful particular cases are: 

Aoo 



Bo(g',0,0) = 
Bo{qlMlM') = 

So(?^0,M2) = 



Bo{qVl^') , 



167r^ 

1 f 1 — 

__ JA^ + ln-2- + l| + 5o(g^M2), 



(A.6) 



167r 



2 ^ Aoo + In 



M' 



+ J(g^0,M2), 
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with the finite parts 



B„(^^ M') = J(q\ M^ M^) = {2 - In (^^) } . (A.7) 



where gm = yJl-AM^/q^. 

The relevant three-propagator integrals are: 



Co(g^O,0,M^) = -_l-,|Li,(^l + ^J-Li,(l)|, 

Co(g^M^M^O) = T^^ln^f^^), (A.8) 

iOTT g V(7m + J-/ 

where 

Li2(y) = - ^ In (1 - xy) = - /' - In (1 - x) (A.9) 
JO a: J X 

is the usual dilogarithmic function. 

Appendix B: Lorentz Structures in the Vector Propagators 

In momentum space, the bare vector-field propagator can be written in the form 

{V^^V<^)o = i^^^''"'{q) = -tJ^ A^^'^'^iq) + n>^^''"^{q) , (B.l) 

My q My 

with the antisymmetric tensors 

\^,pa{q) = ^ [gt^pq^qa - 9pyqt,qa -{p^<y)] , 

^nu,pa{q) = [gupQuqa " Qpvquqa - q^QupOva -{p^(j)\ , (B.2) 

Zq 

^HU,pa = 2 ^dupg^cr ~ 9lJ.crgi>p) i 



2 



obeying the properties: 

Q-A^A-^^Q , A-A^A , = Q , ^ + Q = X, 

q^^nv,pa{q) = q" ^,xu,pa{q) = q''^,xu,pa{q) = q" ^nv,p<j{q) = O- 
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(B.3) 



For any antisymmetric tensor Hfj^i^^ab-, the operator Tcd,af3 acts like the identity, i.e. 

H -I = I ■ H = H . (B.4) 
We can then define the antisymmetric inverse G"'^''^ , which satisfies 

H^.,ab G""-'"^ = G^,,al> H^"'"" = . (B.5) 
The inverse propagator in momentum space is given by 



A^'^'^iq) - i I^'^'f" . (B.6) 



Appendix C: Form Factors with Resonances in the Final State 

We present here a few examples of current matrix elements with external resonances. 
They show that in order to implement a correct short-distance behaviour, one needs to 
introduce additional interactions with more than one resonance field. Moreover, the new 
chiral couplings can be easily determined. 



C.l Axial form factor to Sj^q^ 

Let us consider the two-point correlation function of two axial currents = d'y^^^u, in 
the chiral limit: 

U^iq) ^xjd'x e^^- (0|T J^(O)t) |0) = (-g^^q' + q^q^) UAA{q') • (C.l) 

The associated spectral function Im n^^(i) is a sum of positive contributions corresponding 

to the different intermediate states. At large t, it behaves as a constant. Therefore, since there 
is an infinite number of possible states, the absorptive contribution of a given intermediate 
state should vanish at infinite momentum transfer. 

One can easily check that the minimal RxT Lagrangian of ref. [8], which only contains 
interactions linear in the resonance fields, generates an absorptive S^^qK~ contribution with 
the wrong behaviour at large momenta: lmIlAA{t)\sn ^ constant. The problem can be easily 
identified analysing the corresponding form factor, defined through the matrix element 

( 7r-| JtmTs^^IO ) = - 2 i J^sM^) - ^) pI , (C.2) 
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a) 



b) 



c) 



Figure 8: Tree- level contributions to Ts-K{(f)- 

where = {p.,^ +Ps)^- The lowest-order calculation with the R^T Lagrangian (diagrams |H1 a 
andlHlb) gives a constant form factor, 

^sAq') = ^ , (C.3) 
which obviously is not vanishing at infinite momentum transfer. 

The correct large energy behaviour can be recovered adding the interaction term [41] 

4a = Af^({V^S,V}«^}), (C.4) 
which modifies the form factor (diagram |H|c): 

,c.5) 

Imposing that the form factor must vanish as — » oo, the coupling Af"^ is constrained to 
take the value 

Af^ = -4^ = -^. (C.6) 
' Fa V2 ^ ^ 

The resulting form factor adopts then the usual monopolar form 

^sA<l') = ir^- (C.7) 
ivi^ — q 



C.2 Vector form factor to R^^^ R (R = S, P) 

The two-point correlation function of two vector currents has a similar behaviour at short 
distances. Its spectral function behaves as a constant at large momentum transfer, implying 
that the form factors associated with each intermediate state should vanish at ^ oo. 

The minimal R^T Lagrangian generates the matrix elements: 

{Rti{pi)R-{p2)\dj''u\0) = V2{p2-pirJ'RR{q'), (C.8) 
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a) 



b) 



Figure 9: Tree-level contributions to J-'RR^q^). 



where R = S,P stands for a scalar or a pseudo-scalar resonance. The corresponding form 
factors are just constant at lowest order (diagram El a): 

J'RRiq') = l. (C.9) 

This constant contribution originates in the chiral connection of the resonance kinetic La- 
grangians. 

It is possible again to recover the right QCD short distance behaviour by adding the 
interaction terms 

^VRR = ^A^^^(l^^^V^i?V.i?), (C.IO) 
which change the form factors to (diagram Elb) 

^rf) = l + ^A^^^-#^. (C.ll) 
V 2 My-q 

Imposing a proper high energy behaviour one gets the constraint 

A^^^ = ^ = l, (C.12) 
and a monopolar form for the form factors 
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